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ABSTRACT 

We follow numerically the time evolution of axisymmetric outflows driven magneto-centrifugally from the 
inner portion of accretion disks, from their launching surface to large, observable distances. Special attention 
is paid to the collimation of part of the outflow into a dense, narrow jet around the rotation axis, after a steady 
state has been reached. For parameters typical of T Tauri stars, we define a fiducial "jet" as outlined by the 
contour of constant density at 10 4 cm" 3 . We find that the jet, so defined, appears nearly cylindrical well above 
the disk, in agreement with previous asymptotic analyses. Closer to the equatorial plane, the density contour 
can either bulge outwards or pinch inwards, depending on the conditions at the launching surface, particularly 
the mass flux distribution. We find that even though a dense, jet-like feature is always formed around the 
axis, there is no guarantee that the high-density axial jet would dominate the more tenuous, wide-angle part of 
the wind. Specifically, on the 100AU scale, resolvable by HST and ground-based adaptive optics for nearby 
T Tauri winds, the fraction of the wind mass flux enclosed by the fiducial jet can vary substantially, again 
depending on the launching conditions. We show two examples in which the fraction is ~ 20% and ~ 45%. 
These dependences may provide a way to constrain the conditions at the launching surface, which are poorly 
known at present. 

Subject headings: galaxies: active — ISM: jets and outflows — methods: numerical - MHD — stars: formation 



1. INTRODUCTION 
1.1. The Quest for Large-Scale Wind Structure 

In a previous paper i lKrasnopolskv. Li. & B landfordl [l999l 
hereafter Paper I), we described a modification of the ZEUS 
MHD code and its application to the problem of launching 
a magneto-centrifugal wind from a Keplerian disk. The disk 
was idealized as a boundary condition at the equatorial plane, 
on which open magnetic fields are firmly anchored at all 
radii inside the simulation box. Mass was injected superson- 
ically but sub-Alfvenically onto the open field lines at a pre- 
scribed rate and accelerated magneto-centrifugally along the 
field lines to produce a high speed wi nd. Our approach fol- 
lowed that of lUstvugova et alJ i 19951) and lOuved & Pudritzl 
( 1997). It permitted a decoupling of the dynamics of the wind 
from that of the disk and allowed for a thorough investigation 
of the structure of the wind. Previous numerical studies of 
magneto-centrifugal winds were limited to a relatively small 
region not far from the launching surface. In this paper, we 
seek to extend the wind solutions to large, observable dis- 
tances, and to relate the observable flow quantities on large 
scales to the launching conditions. 

To study the large-scale structure of a disk-driven magneto- 
centrifugal wind numerically, a technical problem must 
first be overcome. This problem is illustrated in Fig. ^ 
(taken from Paper I), where a disk-wind solution repre- 
sentative of those found in existing time-d ependent simu- 
lation s (e.g.. iRo manova et alJl!997t [Ouved & Pudritzl [19971 
lUstvugova et all 19991) is shown. Note that a large fraction of 



the wind remains sub-fast-magnetosonic in the computational 
realm. Information from the edge of the simulation box can 
propagate upstream in a sub-fast-magnetosonic region, creat- 
ing an undesirable coupling between the box boundaries and 
the wind-launching. This coupling, we believe, is at the root 
of the sensitive dependence of wind solutions on the size and 
shap e of t he simulati on box fou nd previously by us and others 
(see 32. ll and lUstvugova et alJl999l for a discussion). In some 
cases, merely changing the shape of the box can destabilize 
a wind and induce chaotic flow fragmentation. This depen- 
dence must be eliminated before one can study the large scale 
structure of the magneto-centrifugal wind with confidence. 

One way to overcome the above difficulty is to restrict the 
wind launching to the inner region of an accretion disk. In 
the specific context of star formation, there are reasons to be- 
lieve that the bulk of the outflow is indeed driven from an 
inner disk. There are two popular magneto-centrifugal mod- 
els for the production of optical jets, su ch as the spectacular 
HH 30 jet _ studie d in detail with the HST jBurrows et al J 19961 
Rav et al. 1996). The models differ on where the wind-driving 
open magnetic fields are anchored. If the open field lines are 
anchored in a narrow region of disk near the corotation radius 
of the stellar magnetosphere, then the wind is called an "X- 
wind" iShu et alJl2000l) . The X-wind serves the fundamental 
purpose of removing angular momentum from the central star 
and keeping it rotating at a rate well below the breakup rate, 
as observed. If, on the other hand, the field lines are anchored 
over a wider region of the disk, then a "disk- wind" is produced 
(eg., Konigl & Pudritz 2000). The disk- wind may be primar- 
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FIG. 1 . — A representative magneto-centrifugal wind launched from a Ke- 
plerian disk (in arbitrary units; taken from Paper I). Shown are magnetic field 
lines (light solid lines), velocity vectors (arrows) and the fast magnetosonic 
surface (solid line of medium thickness). The thickest solid line divides the 
portion of the wind that becomes super fast-magnetosonic inside the simula- 
tion box (above) from the portion that does not (below). 



ily responsible for driving the mass accretion through the por- 
tion of the disk where the wind is launched. Wind launch- 
ing should be easier from the inner part of protostellar disks, 
where the temperature is high enough (T > 10 3 K) that ther- 
mal ionization of alkali metals ensures that the disk material 
and wind-launching magnetic field are well coupled. In the 
cooler outer part, the bulk of disk material is too weakly ion- 
ized to participate in the wind launching, althoug h in the sur- 
face l ayer and at large radii (say ~ 10 2 AU; War dle & Konigll 
119931) protostellar X-rays and Galactic cosmic rays (if not ex- 
cluded by the wind) may provide a level of ionization suffi- 
cient for good coupling. We will not consider in this paper 
the possible outer disk wind, which is probably less relevant 
for the fast jets and winds of hundreds of kms" 1 that we are 
interested i n. S uch a wind can modeled with our simulation 
setup (see 32. 1> if desired since the size of the launching re- 
gion is fully adjustable. 

The large-scale structure of magneto-centrifugal winds is 
a prerequisite for quantitatively modeling a variety of out- 
flow phenomena observed around YSOs, including optical 



jets, neutral atomic winds, and bipolar molecular outflows. It 
is the observational data on these phenonema that will con- 
strain the conditions at the wind launching surface, which 
are poorly known. Particularly tight constraints are expected 
from high spectral as well as high spatial resolution obser- 
vations of forbidden line emission from T Tauri winds. A 
recent example of such observations is that of Ba cciotti et alJ 
(2000). They constructed synthetic two-dimensional images 
of the DG Tau wind based on HST/STIS data in several for- 
bidden lines (including [OI]AA6363,6300, [SH]A6731, and 
[NII]A6583) and H a in several radial velocity intervals, with 
a spatial resolution of order 0.1" (or about 15 AU at the dis- 
tance of DG Tau). To take full advantage of such observa- 
tions, one needs a reliable model that can follow the flow from 
the (small) launching surface to observable distances. As a 
first step towards such a model, we restrict our attention to 
those magneto-centrifugal winds that can reach a steady state, 
even thoug h most observed jets a re time-dependent, possibly 
episodic (Reip urth & Ragalll999l) . We postpone an investi- 
gation of the intrinsically time-dependent winds to a future 
paper. 

Although our intended application is primarily to YSOs, we 
note that magneto-centrifugal winds driven from a limited re- 
gion of accretion disks may occur in other astrophysical sys- 
tems, e.g., cataclysmic variables (CVs) and possibly plane- 
tary nebulae. A related jet launching mech anism, t hrough tor- 
sional Alfven wave, was first discussed in Uchida & Shibata 
and refined in man y subsequent works (see, e.g., 
Kato, Kudo, & Shibata 2002 and references therein). 

1 .2. Uncertain Jet Formation Efficiency 

Asymptotic structure o f steady magneto-cent rifugal winds 
has been analy zed by iHevvaerts & Normanl dl989l) and 
IShu et alJ |l995), among others. The focus of such anal- 
yses was on jet formation, through the hoop-stress associ- 
ated with toroidal magnetic field. While there is a gen- 
eral agreement that a dense, narrow jet can indeed form in 
the a xial region of a magneto-centrifugal wind (see, how- 
ever, Qkamotol II1999H for a different view), it is not clear 
how the jet properties, such as their density, mass and mo- 
mentum fluxes, depend on the launching conditions. Indeed, 
the term "j et" was used by different autho rs to mean differ- 
ent things. Bogovalov & Tsinganos] dl999l) . for example, re- 
fer to the axial region with a more or less uniform density 
and magnet ic field strength and vertical streamlines as the 
jet, whereas fShu et al.l dl995 ) suggested that the observed jet 
is simply the dense part of an underlying anisotropic wind 
above, roughly speaking, some fiducial density that shows 
up more prominently than the more tenuous equatorial part 
in the emission of f orbidden fines such as [OI] and [SII] 
(Shan g. Shu. & G lassgoldl ["l998l) . We adopt the latter mean- 
ing of "jet", which is more observationally oriented. For 
simplicity, the fiducial density at the jet boundary is chosen 
to be hh = 10 4 cm~ 3 , comparabl e to the values derived by 
iBacciotti. Eisloffel. & Ravl (f!999) for the famous HH 30 jet 1 . 
An important question is: what fraction of the mass flux of 
the wind is collimated into the fiducial jet? We will address 
this question using numerically determined, large-scale wind 
solutions. 

The fraction of mass flux residing in the jet measures, in 

1 We note that the density derived bv lBacciotti et alj|1999l) for the HH 30 
jet is not constant along the jet, but decreases from a value ~ 10 5 cirT 3 close 
to the exciting source to a value ~ 10 4 cm~ 3 further out. 
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some sense, the efficiency of jet formation. It can be con- 
strained by observations. Take the HH 30 system as an 
example. For this object, one can eliminate a jet forma- 
tion ef ficiency as low as 10% from the following consider- 
ations. iBacciotti et alJ (^999) showed that the mass flux in 
the main jet is ~ 1.7 x lO _9 M yr 1 (~ 1.8 x lO^Mgyr" 1 
in the counter jet) close to the exciting source. If only 10% 
of the mass flux resides in the jet, the total mass flux of 
the wind from both sides would be ~ 3.5 x lO^M^ yr" 1 . 
If a t hird dShu et alJ[l995h or less (e.g., iFerreira & Pelletierl 
119951) of the material accreted through the disk gets ejected 
in the wind, then the disk mass accretion rate must be ~ 
10 Mgyr -1 or higher. Such a large accretion rate would 
be difficult to accommo date by the very weak IR excess 
from ISO observations (Stapelfeldt & Monet] [I9981) . which 
points to a relatively inactive disk, and by the low disk mass 
(~ 10~ 3 Mrri) estimated fro m scattered light and dust emission 
(Stapelfeldt & Padgett 2001), although this estimate is uncer- 
tain. Increasing the jet mass flux fraction to 50% would lower 
the disk accretion rate to ~ 2 x lO^M^yr" 1 (a more typical 
value for T Tauri disks; Gullbring et alJll998» . if a third of 
the accreted material is ejected. If a smaller fraction of the 
accreted material is ejected, then the required accretion rate 
would be higher. It therefore appears that, at least for HH 30, 
the jet formation efficiency must be fairly high, of order 50% 
or more. It is not clear whether such a high efficiency can be 
achieved naturally in a magneto-centrifugal wind. 

The rest of the paper is organized as follows. In 5j2j we 
present our formulation of the problem of launching magneto- 
centrifugal winds from only inner disk regions. It is followed 
by a discussion of a "standard" wind solution which illustrates 
several generic features of jet formation and large-scale struc- 
ture of steady magneto-centrifugal winds in In §J3] we ex- 
amine the effects of the mass flux distribution at the launching 
surface on the jet/wind properties. We conclude and discuss 
future work in |5] 

2. FORMULATION OF THE INNER DISK-DRIVEN WIND 
PROBLEM 

2.1. Implementation of a Finite Launching Surface 

As mentioned earlier, we limit the launching of magneto- 
centrifugal winds from the inner, relatively hot regions of pro- 
tostellar disks. Numerically, this restriction enables us to ex- 
tend the wind solution to large distances without having to 
worry about the undesirable coupling between the simulation 
box and the launching region, which occurs in the sub fast- 
magnetosonic region in the lower-right corner of the solution 
shown in Fig. ^ The reason for that region of flow to re- 
main sub fast-magnetosonic is simple: the wind coming off 
the outer part of the disk encounters the edge of the simula- 
tion box too soon. In other words, the wind in the region sim- 
ply does not have enough room to get accelerated to the fast 
magnetosonic speed. This situation remains true as long as 
the wind is driven off all of the (equatorial) disk plane (as as- 
sumed in all previous time-dependent disk- wind simulations), 
regardless of the box size. This troublesome feature can be 
eliminated by restricting the wind launching to only the inner 
region of a disk (within a radius denoted by R£), as sketched 
in Fig. 13 To fill all available space above (and below) the 
equatorial plane, we demand that the last field line lie exactly 
on the equatorial plane 2 . Wind plasma sliding along the last 

2 Of course, if the wind is confined externally, either by a flaring disk or 
an ambient medium, the last field line must be modified accordingly. Exter- 




surface 

FIG. 2. — Schematic view of a magneto-centrifugal wind launched from a 
limited, inner disk region. 



(horizontal) field line, anchored at the radius Ri, will become 
super fast-magnetosonic in the computational realm, provided 
that the size of the simulation box is sufficiently large. Once 
the whole fast magnetosonic surface is completely enclosed 
inside the simulation box, the size and shape of the box should 
have minimal effects on the structure of the wind, especially 
near the launching surface, since information cannot propa- 
gate upstream in a super fast-magnetosonic region. In this 
way, we should be able to study the wind structure up to arbi- 
trarily large distances from the source region, limited only by 
computer time. 

There is, h oweve r, an additional consideration. 
lUstvugova et alJ \ 1999ft have shown that artificial de- 
pendence of simulation results on computational box shape 
can appear even in a super-fast-magnetosonic flow, if the 
Mach cones at the boundaries intersect the simulation region. 
The opening angle <p of the Mach cone around a streamline 
is given (for a cold flow) by isxp-ip = 1/(M 2 - 1), where 
Mf is the fast Mach number. As the flow becomes more 
super-fast-magnetosonic, the cone gets narrower, more 
forward-directed, and has less chance of intersecting the 
simulation region. Our simulations, designed to produce 
clearly super-fast flows, fulfill this requirement everywhere 
except for a small region around the axis. As we discuss 
below, the axial region receives a special treatment, and does 
not appear to create any numerical problems. 

2.2. Equations 

For completeness and to define symbols, we list below the 
time-dependent MHD wind equations that we solve numeri- 

nal confinement of magneto-centrifugal winds is an important topic that we 
plan to explore numerically in the future. It is expected to play a more domi- 
nant role in the outflows of embedded sources than in the winds of optically 
revealed T Tauri stars. Note that the dense material inside the disk beyond 
the wind launching region prevents the "last" field lines of both sides of the 
disk from reconnecting. Also, the interface between the wind and disk ma- 
terial outside Rl may become Kelvin-Helmholtz unstable. Such additional 
complications are ignored in this paper. 
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cally: 



^ + V-(pv) = 
at 

d\ 

p— + p(\- V)v = -Vp + pV$ s +j xB/c 

— = Vx(vxB) = Vx£ 
at 



On 
dt 



+ V-(mv) = -/?V-v 
p = ( 1 -\)u 



(1) 

(2) 

(3) 

(4) 
(5) 



where 



p = matter density 
v = velocity flow field 
B = magnetic field 
j = (c/47r)V x B = current density 
£ = v x B = -cE where E is the electric field 
$ g = gravitational potential 
p = thermal pressure 

// = internal energy density (per unit volume) 
7 = adiabatic index 

We have adopted a cylindrical coordinate system (R,(f>,z), 
with the spherical radius given by r = y/R 2 +z 2 . The gravity 
field is smoothed inside a sphere of radius r g , through 



n 2 



•3r/), 



(6) 



where £lo and r ? are fixed parameters, whose product 
vko = ^o>> gives the scale for the rotational speed vk = 

R (9$ g /9/?)__ = QqR on the disk inside r g . Outside r g , 

the gravity is not smoothed, and the potential has the usual 
form 



Q, 2 r 3 



It yields a Keplerian speed of 



vk = v Ko 



R 



1/2 



(7) 



(8) 



on the disk. The parameter Slo is simply the angular speed at 
the inner edge of the (unsoftened) Keplerian disk. We label 
magnetic field lines with the magnetic flux function ip defined 
by dip/dR = RB Z . Except for a few specific differences noted 
here and below, the notations and numerical methods we use 
are the same as those in Paper I, to which we refer for more 
details. 

The governing equations are solved numerically using the 
ZEUS code, subject to a set of appropriate initial and boundary 
condit ions. We use a p arallelized version of the Zeus 3D code 
iClarke. Norman. & Fiedle rl 119941). wh ose main algorithms 
are based on Stone & Norman ( 1992a b) and Hawlev & Stone 
The magnetic field is evolved by the well -known 
constrained transport method (Evans & Hawlev 1988), which 
preserves the condition V ■ B = at all times as long as this 
condition is satisfied everywhere initially. 

2.3. Boundary Conditions 

On the four sides of our 2D simulation box (see Fig. [2}, 
we need to impose boundary conditions. Numerically, these 



conditions are imposed by assigning values to flow quanti- 
ties in a few ghost zones just outside the active computational 
grid. The ghost zone quantities include the density p, internal 
energy u, the three components of velocity v, and the three 
components of the field £ = v x B. 

The conditions on the axis R = and the outer boundaries 
R = ^max and z = Zmax are relatively straightforward: on the 
axis we demand reflection symmetry, whereas on the outer 
boundaries we allow for outflow. These conditions are stan- 
dard options for the Zeus code. The conditions on the equa- 
torial plane z = are more difficult to specify and implement. 
They represent the main technical improvement that distin- 
guishes this work from others. 

For our problem, the z = surface is divided into two dis- 
tinct regions on which different boundary conditions are im- 
posed (see Fig. Inside the radius Rl, which divides these 
two regions, magnetic field lines are anchored and outflow- 
ing materials are injected. The boundary conditions there are 
the same as those used in Paper I. Specifically, the quanti- 
ties p, B, and v. are specified as known functions of the ra- 
dius R at all times, and the two horizontal components of the 
£ field satisfy = and £r = v&B z . To implement these 
boundary conditions for £, the horizontal components of v in 
the ghost zones are chosen so that vr/v z = Br/B z , and = 
VK + B^Vp/Bp = vk+B^v z /B z . In the ghost zones below z = 0, 
we force the components of £ to have values: £^(-z) = -£<p(z), 
£ R (-z) = 2v^{R)B z {R)-£^{z), and £ z (-z) = £ z (z). These con- 
ditions allow the foot point of a wind-launching open field 
line to be pinned firmly on the launching surface and at the 
same time leave the field line free to bend both radially and 
azimuthally in response to the stresses exerted by the wind 
material. 

Our boundary conditions o n the launching surface inside 
Rl are similar to those used bv lUstvugova et af] J 1 9951 Tl 9991) 
and Roman ova et alJ dl997l) . except that they also considered 
the case of subsonic wind injection. In the subsonic case, the 
wind must pass smoothly through a sonic point close to the 
launching surface, which fixes the mass loading rate and re- 
duces the number of the quantities that must be specified at the 
z = surface by one. The sonic point is sensitive to the ther- 
mal and magnetic structure of the disk atmosphere, which is 
not well understood. We will postpone a detailed treatment of 
the son ic transition to a future investigation. lOuved & Pudritzl 
i 19971) appear to have fixed all flow quantities in the ghost 
zones, which is different from our treatment. Fixing all quan- 
tities in the ghost zones risks creating discontinuities in the 
radial and toroidal components of the magnetic field between 
the ghost and active zones, since in the active zones these two 
components cannot be prescribed; they are determined by the 
stresses in the wind. 

Outside the radius Rl, the equatorial z = surface is occu- 
pied by the last flux surface along which the wind material 
emanating from the radius Rl flows. We demand that the ver- 
tical field component B z = on this portion of the z = surface, 
and B z (-z) = -B z (z) in the ghost zones below z = 0. These con- 
ditions on the magnetic field are enforced by the conditions on 
the field £, which are £^(-z) = -£<t>(z), £r(~z) = -£r(z), and 
£ z (-z) = £ z (z). Equation Q then ensures that the initial value 
of B z at z = will not change. For the hydrodynamic quanti- 
ties, we demand that p, u, vr and be continuous across the 
z = surface and assign the value v z = to the z = surface 
and the ghost zones below. 

To complete our discussion of boundary conditions on the 
z = surface, we note the following three additional features. 
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First, to fill out space, the field lines anchored in the wind 
launching region must become more and more horizontal as 
the edge of the region (at R = Rl) is approached. This con- 
dition is satisfied by forcing the functions defining v z (R) and 
B Z (R) for (z = 0, R < R L ) to zero at R = R L , keeping their ra- 
tio finite. Second, to prevent inflow in the region surrounding 
the rotation axis where the magneto-centrifugal mechanism 
fails, we inject a cold, light, fast axial flow capable of escap- 
ing the gravitational pull of the central compact object. Plau- 
sible physical justifications for this fast injection are given in 
Paper I. Finally, during the wind acceleration, there can ap- 
pear temporary negative values of v z in active zones next to 
the launching disk. For this case only, we allow the boundary 
condition at z = to absorb the flow, by setting v z to zero at 
z = 0, and allowing the other hydrodynamic quantities to have 
the same value in the ghost and active zones, but not alter- 
ing the conditions on the £ field. This procedure avoids some 
numerical defects, as explained in Paper I. The need for this 
special treatment disappears once a steady state near the disk 
surface is reached. 

2.4. Initial Magnetic Field Distribution 

Following Ouved & Pudritz (119971) . we choose as the initial 
(purely poloidal) magnetic field distribution above the disk a 
current-free field compatible with the function B Z (R) defined 
at z = 0. Such an initial field has the advantage of perturb- 
ing the ambient medium as little as possible. (This gas will 
eventually be pushed outside the computational domain by the 
fast moving wind launched from the disk.) One can solve the 
force free equation V 2 V> = for the flux function ip through 
separation of variables, with the following set of boundary 
conditions: ip(R = 0,z) = 0, ip(R,z = oo) = 0, tp(R = oo,z) = 0, 
and ip(R,z = 0) = fofS). The function i/jd(R) is prescribed in 
the launching region inside the radius Rl', beyond Rl, it has a 
constant value ipL- A formal solution that satisfies these con- 
ditions is 



Tp(R,z) = R / dke- kz Ji(kR)g(k) 
Jo 



(9) 

g(k) = k I dRJi(kR)ip D (K) (10) 
Jo 

i-Rl 

= k / dRJ l (kR)yj D (R) + J (kR L )^ L (11) 
Jo 

where Jo and Ji are Bessel function of zeroth and first order. 
Doing these integrations at each grid point can be expensive, 
especially in regions where the integrals are slow to converge, 
due to the rapidly oscillatory integrand. A more convenient 
and surer procedure to find the solution consists of first eval- 
uating with sufficient approximation the integrals at the outer 
boundaries R = /? max and z = z m ax, and then using these values 
to solve an elliptic problem with Dirichlet boundary condi- 
tions, together with the known values of ipD at the disk and 
axis. Since the size of the simulation box is much larger than 
the size of the launching region, one can expand the integrals 
in powers of k at distances far from the launching region (i.e., 
r ^> Rl). The first two terms of expansion, which proved suf- 
ficient to calculate our boundary values, are 

^(* )2 ) = ^(l-;) + ^£, (12) 

where C = f* L Rip D (R)dR-ip L RL 2 A We solved the elliptic 
equation by a relaxation method, taking as an initial guess the 
simple distribution ip(R,z) = ^D{r){\—z/r). 



With the governing equations, numerical method, initial 
and boundary conditions explained, we are ready to exam- 
ine in detail a representative wind solution that has reached a 
steady state. 

3. AN ILLUSTRATIVE WIND SOLUTION 

3.1. Simulation Setup 

We first establish a "standard" run against which other sim- 
ulations will be compared. All simulations are carried out 
with dimensionless quantities, but to facilitate comparison 
with high resolution observations that are becoming avail- 
able, especially from the HST, we will present results with 
dimensional units. To fix units, we adopt one solar mass for 
the central object M*, and Ro = 0.1 AU for the inner radius 
of the wind-launching region, which is roughly the corota- 
tion radius of the stellar magnetosphere (inside which t he disk 
is either very s ub-Keplerian or non-existent; Ktinigl 1991; 
LShu etal] fl995). The gravitational softening radius r g is set 
equal to Ro- The central mass and inner radius set a scale 
for the velocity, which is the Keplerian speed at Rq, vk = 
^JGM^/Rq = 94kms _1 . The scales for the magnetic field 
strength, Bo, and the density, po, are related by the require- 
ment that Bq/(4ttpo) = vko 2 - Both of these scales are fixed 
by demanding that the total mass flux of the outflow M (from 
each side of the disk) be Mo = 10~ 8 M@ yr" 1 , a value represen- 
tative for classical T Tauri stars (e.g., Edwards, Rav, & Mundt 
1993). 

For the standard run, we choose an outer radius for the wind 
launching region, Rl = 1 AU, ten times the inner radius Rq. On 
the equatorial plane z = inside Rl, one needs to specify the 
distributions of the density p at the base of the wind, and the 
vertical components of both the injection velocity (v,) and the 
magnetic field strength {B z ). We consider a single-component 
distribution for the magnetic field 



B Z (R) = B,,, 



i+(R/R y 



-et/2 



x f(R), 



(13) 



where the parameter B z max sets the scale for the field strength 
and the exponent et the distribution. They are taken to be 
4 Bo and 1 .5 respectively for the standard run (where Bo is the 
magnetic field unit defined above). The function f(R) is set to 
unity between the origin and a transitional radius Rj, and 



f(R)= 1 + 2 



R z -Ri 



Rl 



-Rl 



Ri -R 



Rl 



-R 2 T 



(14) 



between Rj and the outer radius of the launching region Rl- 
It is used to force both B z and v z to vanish together at the 
outer edge of the la unch ing region, as required by the bound- 
ary conditions (see 32. 3> . We choose R j = 0.8 AU. 

For the injection speed at the launching surface, we con- 
sider a two-component distribution, with 



v z (R) = 



l + (R/R ) 



in the fast injection axial region inside Rq, and 
v z (R) = a v K x f(R) 



(15) 



(16) 



in the magneto-centrifugal wind launching region between Rq 
and Rl- The parameter v z o sets the maximum injection speed 
on the axis, and is taken to be 160kms _1 for the standard run. 
The parameter a denotes the ratio of the injection speed to 
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the Keplerian speed, and is set to a small value of 0. 1 for the 
standard run. 

For the density distribution at the base of the wind z = 0, 
we adopt a small, constant value of p m = 0.1 po (where po is 
the density unit) in the fast injection region inside the radius 
Rq. On the Keplerian disk between Rq and Ri, we prescribe a 
power-law distribution 

P(R) = Pantry , (17) 

up to the outer edge of the launching region Ri. The parame- 
ters p out and e p set the density scale and its decline rate with 
radius. They are taken to be p out = 10po an d e p =\ for the 
standard run. Note the factor of 10 2 contrast between the den- 
sities of the fast axial injection (p m ) and the inner edge of the 
magneto-centrifugal wind (p out ). 

The launching conditions specified above yield the follow- 
ing units for the field strength and density: Bq= 1 . 1 G and 
po = 1.2 x 10~ 15 gcm" 3 . These conditions are displayed in 
Fig. together with a distribution of the accumulative mass 
flux injected into the wind and the magnetic flux enclosed 
within a given radius on the disk. It is clear that the fraction 
of mass flux in the fast injection region inside R = Rq is small 
(~ 1.5%), despite the large injection speed in this region; the 
large speed is more than offset by the low density specified. 
The large velocity shear near R = Rq does not appear to cre- 
ate any instability, possibly because of the strong (poloidal) 
magnetic field present in the region. 

Above the equatorial surface z = 0, we start with a purely 
poloidal, potential mag netic field, constructed using the 
method outlined in 32.41 The computational domain is ini- 
tially filled with a low density medium, whose density dis- 
tribution is prescribed either analytically or taken from other, 
similar runs that have reached a steady state. Judicious choice 
of the initial ambient density distribution can speed up the 
flow convergence to the steady state enormously, although the 
final wind solution does not depend on the choice, since the 
ambient medium will eventually be swept out of the simula- 
tion box completely. To isolate the magneto-centrifugal wind 
acceleration and collimation, the focus of our study, from 
thermal effects, we keep the flow cold at all times by arbitrar- 
ily setting a small value of 0.02kms -1 for the sound speed. 

Our main goal is to study the structure of magneto- 
centrifugal winds on large, observable scales. For this pur- 
pose, we adopt a simulation box of 100 AU x 100 AU. It is 100 
times the size of the wind launching region and 1000 times 
the inner Keplerian disk radius Rq. Such a scale is directly 
accessible to high resolution obs ervations using the H ST 
and ground-based adaptive optics (Dougados et al. 2000) for 
nearby optically revealed, T Tauri winds. Outflows from 
deeply embedded sources can be probed on a similar scale us- 
ing water masers. For example, the water masers of the Class 
source SI 06 FIR appear to lie on a 3AU (width) x 4AU 
(length) U-shaped surface some 25 AU away from the cen- 
tral source, possib ly driven by a well-collimated "micro jet" 
iFuruva et alJl99 9). Numerically, we cover the active compu- 
tation domain with 190x210 zones. In the vertical direction, 
there are 40 uniform grid points covering the region from to 
0.8 AU, and 150 non-uniform zones covering the region from 
0.8 to 100 AU, with the zone size increasing by a small factor 
of 1.035023 between adjacent zones when moving outwards. 
In the horizontal direction, the first 60 grid points from to 
1.2AU (thus covering the entire launching region) are uni- 
form, followed by 150 non-uniform grid points with approx- 



imately the same ratio between zones as before (1.034987). 
Larger boxes are in principle possible, but limited by com- 
puter time. 

3.2. Global Structure, Flow Collimation, and Jet Formation 

We follow the evolution of the material injected into the 
computation domain numerically with our modified Zeus 
code, keeping the boundary conditions shown in Fig.|3]fixed 
at all times. The injected wind material pushes aside the low- 
density ambient medium and propagates towards the outer 
edges of the simulation box. After all of the ambient medium 
is swept out of the box, the wind settles quickly into a steady 
state, with no appreciable time variations. The overall appear- 
ance of the steady-state wind is given in Fig.@] on two scales. 
In Fig.|4k, we show on the 10AU scale a meridian view of 
nine streamlines (also field lines) that divide the wind into 
ten zones of equal mass flux. Superposed on the streamlines 
are the isodensity contours (in shades), with values of density 
decreasing exponentially outwards. A salient feature is the 
fast magnetosonic surface, shown in dashed lines in Fig. 0^. 
As mentioned earlier, the fast surface closes inside the sim- 
ulation box, except near the axis where the light, fast axial 
flow resides. The tenuous flow in this narrow axial region 
remains dominated by the (largely poloidal) magnetic field, 
which can, in principle, provide some stability to the wind as 
a whole against kink instabilities (Shu et al. 1995). The axial 
injection occupies a decreasing fraction of the total volume 
on larger scales. It becomes nearly invisible in panel (b) of 
Fig.|4] where streamlines and isodensity contours are plotted 
on the 100AU scale. Here, most of the space is filled with 
the flow driven from the disk magneto-centrifugally, whose 
acceleration and collimation properties we want to study in 
detail. 

We first concentrate on the collimation properties of the 
standard solution. Self-collimation is commonly accepted as 
a hallmark of the magneto-centrifugal winds. The degree 
of collimation i s still a matter of debate, however (see, e.g., 
Okamoto 1999). Our self-consistently determined solution 
up to large distances can shed light on this important issue. 
From panel (a) of Fig. |4] we find that the standard wind is 
not well collimated on the small scales close to the launch- 
ing surface, although hints of gradual collimation are evi- 
dent in the shapes of both streamlines and density contours: 
the streamlines bend slightly towards the rotation axis, which 
forces the density contours to elongate along the axis. The 
elongation of density contours shows up more clearly on the 
larger scale in panel (b) of Fig. |4] with the contours appear- 
ing nearly cylindrical in the region within ~ 15° of the ro- 
tation axis. Asy mptotically cylin drical density stratification 
was predicted bv lShu et alJ (fT9 95) analytically. The cylindri- 
cal shape does not extend all the way to the equator, however, 
at least on scales up to 100AU. The density contours in the 
equatorial part of the wind bulge o utwards. The b ulging is 
similar to mat found nu merically by S akurail d!987l see also 
lOuved & Pudritzlll997l) . The relative shape of the bulge ap- 
pears to change little from one isodensity contour to the next, 
which is probably related to the logarithmically slow collima- 
tion ofstreamIinesJnJhe_£gra region at large distances 
(see Bogovalov & Tsinganos 1999 for a recent discussion). It 
will be slow to disappear, if ever, on even larger scales. 

The shape of isodensity contours on large scales is impor- 
tant to compute because it is intimately related to the ob- 
served shape of optical jets. A density stratification bulging 
out prominently near the base as in Fig.|4]may be in conflict 
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Cylindrical radius (AU) Cylindrical radius (AU) 

FIG. 3. — The conditions specified on the wind launching surface z = for the standard run in dimensionless units. The dimensional units are vko = 94kms~' 
for velocity, Mo = 10~ 8 Mq yr~' for mass flux, Bo = 11 G for field strength, and po = 1.2 X 10~ 15 gem -3 for density. The magnetic flux ip has a unit of BqRq 2 , 
where Rq = 0. 1 AU is the length unit. 





FIG. 4. — Streamlines (light solid) and isodensity contours (heavy solid lines and shades) of the standard steady wind solution on the (a) 10 AU and (b) 10 2 AU 
scale. The streamlines divide the wind into ten zones of equal mass flux. The dashed line is the fast magnetosonic surface. The arrows are for poloidal velocity 
vectors, with length proportional to the speed. There are two shades per decade in density, and the densities labeled are the number densities of hydrogen nucleus 
for a wind of 10~ 8 Mq yr~' (per hemisphere), assuming a helium abundance of 10% in number. 



with high resolution observations of some H H jets. In par- 
ticular, the spectacula r HST images of HH 30 tB urrows et al.l 
ll996tlRav et all 19961) show that this jet remains nearly cylin- 
drical all the way to the disk surface; if anything, it appears to 
be slightly pinched in the equatorial region, keeping in mind, 
however, that contamination of the jet emission from scattered 
stellar light is a concern near the disk. Since the emissivity is 
sensitive to density, the optical appearance of a jet should fol- 
lowto a large extent that of isodensity contours ( Sha ng et al.l 
1998), although the thermal structure affecting the emission is 
uncertain. The requirement that the density contours be nearly 
cylindrical all the way to the disk surface (assuming the likely 



situation of a much less rapid variation in temperature than 
in density) imposes a constraint on possible combinations of 
magnetic field and mass flux distributions on the launching 
surface. The standard run does not appear to satisfy this con- 
straint. 

Another generic feature predicted by the asymptotic anal- 
ysis of Sh u et alJ (1199 5) is the strong stratification in density 
transverse to the rotation axis. This feature shows up clearly 
in the density contours. We stress that the stratification is 
a most important characteristics of the magneto-centrifugal 
wind, even though individual streamlines appear hardly bent 
at all. We estimate that, at a height of 90 AU, the density dis- 
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tribution decreases rapidly outside a small core of ~ 2AU, 
roughly as a power-law of the cylindrical radius ~ R~ l 3 . The 
decline is slowe r than the asymptotic estimate of ~ R~ 2 of 
LShu et alJfl995h . obtained assuming a constant flow speed on 
all field lines; the difference is partly accounted for by the 
lower speed achieved along the field li nes in the more equa - 
torial part of the standard solution (Matzner & McKee 1999). 
Nevertheless, the fact that the density decline is significantly 
faster than R' 1 means that the central high-density "jet" will 
stand out in maps of integrated emission along the line of 
sight. The steep density stratification, coupled with a mod- 
erate anisotropy in flow velocity, gives rise to a strong con- 
centration of energy and momentum fluxes in the direction 
along the axis, with implications on th e dynamics of bipo- 
lar molecular outflows (Li & Shu 1996), which are probably 
formed when the highly anisotropic winds run into the am- 
bient medium. The decrease in density away from the axis is 
qualitatively c onsistent wit h the HST observations of the wind 
of DG Tau bv lBacciotti et all ( 120001) . 

Even though a dense, nearly cylindrical "jet" is naturally 
produced around the axis, at least well above the disk, the 
fraction of the wind mass flux that gets collimated into the 
jet is a concern. To fix ideas, we adopt a fiducial number 
density of «h = 10 4 cm" 3 to delineate the outer boundary of 
the central "jet". From Fig. |4j), we find that only ~ 20% of 
the wind mass flux resides in the fiducial jet at the height of 
z = 10 2 AU; the majority remains in the wide-angle compo- 
nent. Such a jet formation efficiency may be too low, at least 
for the famous HH 30 system, as discussed in the introduc- 
tion. This example illustrates that there is no guarantee that 
a disk-driven magneto-centrifugal wind would automatically 
(a) concentrate most of its mass flux into an axial jet or (b) 
produce a density stratification cylindrical all the way to the 
equator. 

The relatively low efficiency of jet formation and undesir- 
able shape of the density contours of the standard solution are 
related to the particular distribution of mass flux chosen. It is 
dominated by the mass injected at the outer part of the launch- 
ing region, since the mass flux per octave in radius is 



between the inner radius of the Keplerian disk Rq and the 
radius Rj (beyond which the injection speed is modified by 
equation 1141 ). Mass from the outer part is loaded onto more 
horizontally inclined field lines, which are harder to collimate. 
The above deficiencies are rectified to some extent in the wind 
solutions to be discussed in the next section. Moderately col- 
limated winds like the standard solution may be applicable 
to the outflows from young high-mass stars, which could be 
more dominated by t he wide-an gle component than their low- 
mass counterparts ( Rich er et al J20001) . 

We note that the standard solution (and other solutions to 
be discussed below) does not recollimate towards the rota- 
tion axis at large distances, and does not show any oscil- 
lations along the axis. Both behavio r s are seen in some 
self-similar (ICha n & Henriksenl Il 980t iBlandford & Pavnd 
119821 lAchtenbe rg. Blandford. & Goldreichll983l) and width- 
averaged jet solutions (e.g., Spruit, Foglizzo, & Stehle 1997). 
They are most likely a result of the self-similarity assump- 
tion or the averaging over the jet cross-section, not the general 
properties of magneto-centrifugal winds. 

3.3. Flow Acceleration and Kinetic-Energy Domination 



We now concentrate on the acceleration, rather than colli- 
mation, of the standard solution. Magneto-centrifugal accel- 
eration along field lines is illustrated in panel (a) of Fig. [5] 
where the poloidal velocity is plotted as a function of (log) 
spherical radius along the four streamlines that divide the 
wind into five zones of equal mass flux. As expected, the 
flow along each field line crosses the fast point smoothly, get- 
ting accelerated to a speed ~ 2 - 3 times the Keplerian speed 
at the foot point of the field line. Note that for this particular 
solution, the flow moves faster along the streamlines in the 
polar region than in the equatorial region, by a maximum fac- 
tor of ~ 4. This pole-equator velocity contrast is qualitatively 
similar to t hat deduced for the w ind of DG Tau from HST ob- 
servations (Bacciotti et al. 2000). These authors inferred that 
the high velocity component of forbidden line emission of DG 
Tau wind occurs near the axis, and is spatially bracketed by 
the lower velocity component. Whether wind solutions like 
our standard run can explain the density distribution and ve- 
locity field observed quantitatively remains to be seen. 

At large distances from the launching region, most of the 
energy extracted magnetically from the disk, as measured by 
the Poynting flux, is expected to be converted into the kinetic 
energy of the ordered plasma flow. This kinetic domination 
is indeed seen in the standard solution, and is illustrated in 
Fig-EJ>, where contours of constant fast (magnetosonic Mach) 
number, M f = v p /v f (where v p and Vf = [(B 2 p + B 1 ^) / (Air p)] x I 1 
are the poloidal and fast magnetosonic speed) are plotted. 
Note that the fast number along each field line reaches a 
maximum value between ^1.6 and ~ 3.0 at the 10 2 AU dis- 
tance from the origin. Since the ratio of Poynting and ki- 
netic energy fluxes is approximately 2 /Mf 2 at large distances 
(Spruid ll99M) . some — 60-80% of the total energy is ki- 
netic. Similarly high (low) kinetic (magne tic) energy frac- 
tion has been found numeri cally bv lOuved & PudritzM 19971) 
and Ustvu gova et alJ ( 11999ft . This efficient conversion is of 
interest because the remaining, low magnetic energy makes 
potential destruction of tor oidal m agnetic field by kink insta- 
bilities (Spruit 1996; Begelman 1998) less disruptive to the 
wind; the wind has become basically ballistic by the time the 
outer edges of the simulation box are reached, particularly in 
the region not far from the polar axis where the fast number 
is the highest. The transition of the flow from magnetic to 
kinetic domination is an esthetically pleasing feature of the 
magneto-centrifugal wind. The asymptotic kinetic domina- 
tion is in contrast with that of the "classical" 1 -D Weber-Davis 
( 1967) wind solution, where th e kinetic energy flux is at most 
a third of the total dSpruitl 19961) . It also differs from the ultra- 
relativistic case where kinetic domination is difficult, if not 
impossible, to achieve in an idea l, steady MHD wind that 
fills the entire 4 ir steradians dChiueh. Li. & B egelmaiJ H998t 
iQEmoto 2001. 

We should emphasize that even though the fast number Mf 
is greater than unity, it is not greater by much; it has values 
of only a few on the 100 AU scale, as seen from Fig. |5j3. In 
other words, the magnetic field, even though its energy frac- 
tion is small, is not completely negligible on large, observable 
scales. In particular, the relatively modest values of Mf make 
the effects of magnetic cushion important in the interaction 
of the magnetized wind with its ambient medium (or with it- 
self if the wind is variable) through shocks, especially if the 
shocks are oblique. 

4. JETS WITH BETTER COLLIMATION AND HIGHER 
FORMATION EFFICIENCY 
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FIG. 5. — Wind acceleration from small to large distances. Plotted are (a) the poloidal speeds as a function of spherical radius along four field lines that divide 
the wind into five zones of equal mass flux, and (b) the contours of constant fast Mach number Mf , showing that the wind becomes kinetically dominated at large 
distances. Also shown in panel (a) are the fast surface (dashed line) and the location (cross) along each field line where the poloidal speed matches the Keplerian 
speed at the foot point. 



In the standard wind solution discussed in the last section, 
only a relatively small fraction of the wind mass flux resides 
in the fiducial "jet", and the isodensity contours bulge out in 
the equatorial region. Both features are undesirable for mod- 
eling jets like HH 30. In this section, we seek to rectify the 
situation, by concentrating more mass loading at smaller disk 
radii, through a steeper power-law distribution for the density 
at the launching surface while keeping the magnetic field dis- 
tribution fixed. The density distribution is controlled by the 
exponent e p in the equation fl!7t . The standard run has a rel- 
atively flat density distribution specified by e p = 1 . We have 
varied this exponent over a range of values, and found that the 
collimation properties of the wind can change substantially 
from the standard run. We choose a case with e p = 3 to illus- 
trate the changes. 

For the e p = 3 run, we keep the distributions of the mag- 
netic field and injection speed on the disk the same as in the 
standard run, and set the density at the inner edge of the Ke- 
plerian disk p out = 10po as before. Since the density now de- 
clines more rapidly with radius, the (dimensionless) mass flux 
ejected from the Keplerian disk is reduced. We correspond- 
ingly reduce the density p m at the base of the fast injection 
region by a factor of 4 (to p m = 0.025 po), so that its mass flux 
fraction remains small (~ 2% of the total). The vast majority 
of the wind material is launched magneto-centrifugally from 
the surface of the Keplerian disk, according to the distribution 



dM 
dlnR 



-3/2 



(19) 



Nearly half of the mass flux is concentrated at the inner edge 
of the Keplerian disk, between IRq and 1.5 Ro (or 0.10 and 
0.15AU). The concentration is shown clearly in panel (a) of 
Fig. |6] where nine streamlines (also field lines) dividing the 
wind into ten zones of equal mass flux are plotted, after a 
steady state is reached. Most of the streamlines emanate from 
close to the inner Keplerian disk, creating an impression of 
an "X-wind", although most of the wind-launching magnetic 
flux lies outside the region where the mass flux is concentrated 



(see Fig.[3ji). 

The density contours in the e p = 3 case started out more col- 
limated in the launching region than those in the standard run. 
On the 10 2 AU scale shown in panel (b) of Fig. [6] they become 
nearly cylindrical all the way to the equator; if anything, the 
contours pinch slightly inwards in the equatorial region, in 
contrast with the bulging observed in the standard solution. 
Thus, a better collimated jet is produced. 

The efficiency of jet formation is also improved. From 
Fig-IU>, we estimate that ~ 45% of the wind mass flux is en- 
closed, at the height of 100 AU, within the density contour of 
10 4 cm" 3 , which marks the boundary of the fiducial jet. In 
contrast, the jet fraction is only ~ 20% for the standard run. 
Therefore, compared with the standard solution, the wind is 
more dominated by the axial jet, although there is still more 
mass flux in the wide-angle component outside the jet than 
inside the jet itself. For the chosen size of the Keplerian disk 
(one decade in radius, from 0.1 to 1.0 AU) and magnetic field 
distribution (B z roughly proportional to R~ 3 ^ 2 ), we find it diffi- 
cult to collimate more than half of the mass flux into the jet, at 
least on the 10 2 AU scale. The jet fraction may be increased 
by external collimation, which can act in addition to the in- 
trinsic self-collimation. This will be particularly true during 
the earliest, C lass stage of star formation, when a m as- 
sive envelope dAndre. Ward-Thompson. & B arson vl UoOO) is 
present, which can confine the wind. The combination of 
external and intrinsic collimation may produce a more jet- 
dominated magneto-centrifugal wind; such a wind may be 
needed to explain the kinematics of some bipolar molecular 
outflows jLee et alJ2 000). We will address the issue of exter- 
nal collimation in the future. 

We note that the cylindrically-stratified density decreases 
rapidly away from the rotational axis, roughly as R~ 2 2 at a 
height of 90 AU in this example. The decline is steeper than in 
the standard case, where p oc R~ 13 roughly. It is not far from 
R~ 2 , perhaps signaling that the wi nd has rea ched an asymp- 
totic state close to that analyzed by Shu et al. ( 1995). We note 
that the asymptotic wind speed in this particular example is al- 
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FIG. 6. — Streamlines (light solid) and isodensity contours (heavy solid lines and shades) of a steady wind solution on the (a) 10 AU and (b) 10 2 AU scale, with 
a mass loading more concentrated near the inner edge of the Keplerian disk than the standard case. The streamlines divide the wind into ten zones of equal mass 
flux. The dashed line is the fast magnetosonic surface. The arrows are for poloidal velocity vectors, with length proportional to the speed. There are two shades 
per decade in density, and the densities labeled are the number densities of hydrogen nucleus for a wind of 1CT 8 Mq yr~' (per hemisphere). 



most independent of the polar angle, showing little signature 
of the Keplerian rotation at the launching surface. It illustrates 
the obvious point that a disk-wind does not have to be slower 
in the equatorial region; it depends on the distribution of mass 
loading for a given magnetic field configuration. 

As a final remark, we stress that the mass loading in the ax- 
ial region where the magneto-centrifugal mechanism fails is 
somewhat artificial. Nevertheless, it provides a relatively non- 
intrusive inner boundary to the magneto-centrifugally driven 
outflow. It has the added advantage of possibly representing 
a fast stellar coronal wind emanating along open stellar field 
lines; such a wind may have a predominantly poloidal mag- 
netic field, which could add to the stability of the entire flow 
IShu et al.lll995l) . We have done several tests to ensure the 
non-intrusive nature of the fast injection, including increasing 
by a factor of 4 the density in the injection region for the two 
examples discussed in this and last sections. The flow struc- 
ture in the magneto-centrifugal region outside the axial injec- 
tion remains little affected, especially at large distances. Even 
though the axial injection appears to have minimal effects on 
the region of interest, it controls the timestep of the simulation 
through the Courant condition, because of its low density and 
large magnetic field strength (and thus high Alfven speed). 
To explore parameter space in a reasonable amount of time, 
as we do in a companion paper (J. Anderson et al. 2003, in 
preparation), one can in principle adopt a relatively heavy ax- 
ial injection. 

5. SUMMARY AND FUTURE WORK 

We have developed a ZEUS-based numerical code capa- 
ble of following the acceleration and collimation of magneto- 
centrifugal winds from their launching surface to large, ob- 
servable distances. The code is made possible by treating the 
wind-driving accretion disk as a boundary, which allows for a 
clear separation of the dynamics of the disk and the wind. By 
limiting the launching region to the inner part of the disk, we 
were able to follow the outflow up to, and well beyond, the 



fast magnetosonic surface for most of the streamlines. Close 
to the rotation axis where the magneto-centrifugal mecha- 
nism fails, a light, fast jet is injected, which could stay sub- 
fast magnetosonic in the computational domain; it occupies 
an increasingly small fraction of space at large distances, 
and provides a non-intrusive inner boundary to the super-fast 
magneto-centrifugal wind that we are interested in. 

We find that, as expected, most of the magnetic energy ex- 
tracted from the disk can be converted into the bulk flow ki- 
netic energy at large distances, producing an essentially ballis- 
tic wind. Its fast magnetosonic Mach number remains mod- 
erate, however, which has implications for wind interaction. 
Our unique simulation setup allows us to investigate quanti- 
tatively the large-scale wind structure, particularly jet forma- 
tion. I n agreement with the asymptotic analysis of Sh u et alJ 
(1995), we find nearly cylindrical stratification in the wind 
density close to the axis at relatively large heights. Closer to 
the equator, the density contours can either bulge outwards or 
pinch inwards, depending on the launching conditions. The 
fraction of the wind mass flux residing in the dense axial 
jet can vary substantially, again depending on the conditions 
on the disk, particularly the distribution of mass loading rate 
for a given magnetic field configuration. It appears possible 
to have, at least on the 10 2 AU scales around T Tauri stars, 
both winds dominated by the wide-angle component and more 
jet-dominated winds. The exact demarcation of these two 
regimes on the 10 2 AU and larger scales is unclear. It will re- 
quire a time-consuming, systematic exploration of parameter 
space. 

Our robust determination of the large-scale structure of 
magneto-centrifugal jets and winds opens up several new pos- 
sibilities. First, it avoids the artificial instabilities seen in 
some of previous simulations where a major fraction of the 
wind leaves the simulation box with a sub-fast speed. How- 
ever, non-axisymmetric instabilities of physical origin, such 
as the kink instability, are possible. They will be examined 
in a subsequent paper in this series with 3D calculations (R. 



MAGNETO-CENTRIFUGAL LAUNCHING OF IETS: II 



11 



Krasnopolsky et al. 2003; Paper III) using axisymmetric solu- 
tions as a starting point. Second, the size of our wind launch- 
ing region is adjustabl e. It allows us to study both the nom- 
inal "d isk-wind" (e.g. iKonigl & Pudritzl 120001) and the "X- 
wind" (Sh u et aLll2000h . which differ mainly in the width of 
the launching region. The effects of the size of the launching 
region and other parameters, such as the mass loading rate, 
on the large-scale wind structure will be explored in a com- 
panion paper (J. Anderson et al. 2003, in preparation). Finally, 



by adopting a time-dependent launching condition, we will be 
able to study the formation of internal shocks, which may be 
used to interpret the knots commonly observed in YSO jets. 
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